Anomalous roughening of wood fractured surfaces 
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Scaling properties of wood fractured surfaces are obtained from samples of three different sizes. 
Two different woods are studied: Norway spruce and Maritime pine. Fracture surfaces are shown 
to display an anomalous dynamic scaling of the crack roughness. This anomalous scaling behavior 
involves the existence of two different and independent roughness exponents. We determine the local 
roughness exponents (i oc to be 0.87 for spruce and 0.88 for pine. These results are consistent with 
the conjecture of a universal local roughness exponent. The global roughness exponent is different 
for both woods, £ = 1.60 for spruce and £ = 1.35 for pine. We argue that the global roughness 
exponent £ is a good index for material characterization. 



I. INTRODUCTION 

Since the pioneering work of Mandelbrot et. al. |Q, 
it has been firmly established that topography of frac- 
ture surfaces exhibits remarkable scaling properties. A 
fracture surface z(x, y) is statistically invariant under an 
anisotropic scaling transformation: 



(x,y,z) -> (Xx, Ay, A c z) 



(1) 



where £ is the roughness exponent. Experimental results 
obtained on various materials (steels B, glass rocks 
HU, ceramics ]|,|], metallic alloys |6f^|] and aluminium 
alloys |9|,10|]), both fragile and ductile, have shown that 
the roughness exponent £ is found between 0.7 and 0.9 
(sec jllj for a recent review). The robustness of the re- 
sults seems to support the idea suggested by Bouchaud 
et. al. jlo| that £ ~ 0.8 might be a universal value of 
the roughness exponent, i.e., independent of the mate- 
rial properties. This conjecture implies that the fracture 
toughness is not correlated to the roughness exponent 
£. However, morphology of fracture surfaces seems to be 
affected by material properties. 

It has been suggested by Bouchaud et. al. that 
models of front lines propagating through randomly dis- 
tributed impurities fllS-ITq] might be relevant to under- 
stand the morphology of the fracture surfaces p6| , |T7[ . 
The development of the fracture roughness has been de- 
scribed as a Family- Vicsek scalingjlj,^]. However, in a 
very recent experimental study p9|, it has been found 
that the surface of a brittle fracture in a granite block 
exhibited anomalous dynamic scaling properties akin to 
what occurs in some models of nonequilibrium kinetic 
roughening ^-^|. 

The anomalous scaling is defined as follows. The devel- 
opment of the fluctuations of the height h{x, t) with time 
is characterized by the root mean square w(l,t) at time 
t over a window size I along the x-axis (perpendicular to 
the propagation direction) 



w(l,t) 



1 1 1 1 



1/2 



where the brackets denotes an average over the win- 
dow position j. The roughness w(l, t ) is expected to scale 
in the case of anomalous scaling as |24fl 



w(l,t) 



if I < t 1 ! 7 - 
if I > i 1 / 2 



(3) 



where the exponent /?» = (£ — Cloc)/z is an anomalous 
time exponent. This anomalous dynamic scaling involves 
two different and independent roughness exponents: the 
local roughness exponent £j oc , which describes the scaling 
when one considers windows smaller than the system size, 
and the global exponent C for scaling involving the system 
size |24|] . The local roughness exponent Q oc is actually at 
reach of the methods currently used for experiment anal- 
yses. The global exponent £ is more difficult to extract 
from a classical roughness measurement. Both exponents 
have to be taken into account for a complete description 
of the scaling behavior of the surface. According to Eq. 
(||), the correlation length ~ t x l z corresponds to a 
characteristic lenght below which the surface appears as 
self-affine with the local exponent Q oc . 



(2) 



In Ref. |19|, global and local roughness exponents, 
( = 1.2 and Q oc = 0.79 respectively, have been mea- 
sured . The latter study was performed on a me- 
chanically isotropic material (granite). However, many 
materials have anisotropic mechanical properties like 
wood, reinforced concrete and most composite materials. 
Anisotropic properties result generally from structural 
reinforcements along specific directions. It is of great in- 
terest to understand how fractures in such materials are 
influenced by the anisotropic texture. 

In this study, we determine the complete scaling behav- 
ior of the fracture roughness resulting from stable crack 
propagation in wood samples of different sizes. For two 
different woods (Maritime pine and Norway spruce), we 
show that the local fluctuations of crack surfaces exhibit 
anomalous dynamic scaling properties. The global rough- 
ness exponent is different for both woods. Local rough- 
ness exponents are identical for both woods and support 
the conjecture of a universal local roughness exponent 
for brittle fracture surfaces. The main consequence of 
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this anomalous scaling is that the magnitude of the sur- 
face fluctuations over regions is not just a function of the 
region size but also of the system size. 

The paper is organized as follows. In Sec. || we 
describe experimental setups for crack propagation and 
fracture surface measurement. Section fill is devoted to 



the anomalous dynamic scaling behavior. In Sec. IV 



we 



study the roughness magnitude as function of the system 
size. Finally, we discuss implications for fracture process 
in Sec. [v| 



II. EXPERIMENT 

Wood is a natural material which displays a structural 
anisotropy resulting from the presence of running cells 
in the radial direction. Two commercially wood species 
have been tested : Maritime pine (Pinus pinaster Ait) 
and Norway spruce (Picea abies L.). Pine specimens 
have an average oven dry specific weight (p) of 560 kg/m 3 
and growth rings are approximately 4 mm wide. Typical 
values for spruce specimens are: (p) = 390 kg/m 3 and 
growth rings of 2 to 5 mm wide. Moisture content of all 
specimens was measured between 11 and 13 %. 

Crack surfaces are obtained from a modified Tapered 
Double Cantilever Beam (TDCB) specimens. A fracture 
was initiated from a straight notch machined with a band 
saw (thickness 2 mm) and prolonged on few millimeters 
with a razor blade (thickness 0.2 mm). Fracture is ob- 
tained through uniaxial tension with a constant opening 
rate (Fig.0) . The tapered shape of the specimens allows 
to obtain a mode I stable crack growth (see Ref. |25| for 
details) which induces a constant crack speed. The crack 
speed was around 0.6 mm/s (from 0.3 mm/s for small 
specimens to lmm/s for large specimens). Crack sur- 
faces were generated along an average radial-longitudinal 
plane by aligning the growth rings perpendicular to the 
straight notch. In order to obtain an evolution of the am- 
plitude of the roughness as function of the system size, 
three geometrically similar specimens of sufficiently dif- 
ferent sizes have been fractured. We used samples of size 
L equal to 11.25, 30, and 60 mm (see Fig.|l|). 

Anatomical characteristics of wood introduce typical 
scales which might appears as cutoffs for scale imparl- 
ances. Most tetragonal tracheid cells in pine and spruce 
are about 25 pm wide. During loading cell walls break 
revealing U-shaped profiles with rugged edges because of 
the rectangular shape of the tracheid section. Thickness 
of cell walls varies from 2 to 10 pm. This facies of fracture 
surface is characteristic of a local brittle fracture process. 

Topographies of the crack surfaces were recorded with 
a mechanical profiler along regular grids. Grid axes are 
along the x direction which is parallel to the initial notch 
and along the y direction which is the crack propagation 
direction (Fig.pl) . The step of sampling in the x direction 
is adjusted to the minimum cell width : Ax = 25 pm 
and to the cell length in the y direction : Ay = 2.5 mm. 



Profiles along the x-axis were sampled with 2050 points 
for specimens of width L = 60 mm, 1030 points for speci- 
mens of width L — 30 mm, and 360 points for specimens 
of width L — 11.25 mm. For each map, the first pro- 
file (y — 0) is sampled in the immediate vicinity of the 
initial straight notch and has a zero roughness. As the 
distance y to the notch increases, the roughness develops 
up to 3 mm. The vertical resolution is estimated from the 
height differences between two successive sampling along 
the same line. Its magnitude is about 3 pm. Horizontal 
resolutions along x and y-axis are about 5 pm. In the 
case of pine an additional specimen size was tested : L 
= 22.50 mm with 800 points, but only profiles far from 
the notch have been recorded. Table | lists parameters 
of the various studied samples. 



III. ANOMALOUS DYNAMIC SCALING 



As mentionned above, fractures of all specimens have 
been obtained at constant crack speed. Subsequently, we 
assumed a linear relationship between the y-position of 
the profiles and the crack propagation time t. Height 
profiles are considered as descriptions of the advancing 
crack front h(x,t). Complete spatio-temporal evolution 
of the crack front can thus be produced from roughness 
maps. 

In Figure ^J, we present the development of the rough- 
ness w(l,t) versus time t in a log-log plot for different 
window sizes I in the case of the spruce specimen (s60- 
1) which is 60 mm wide. The upper line is a fit of the 
roughness growth for the largest window size (I = 13.975 
mm). The slope of this fit provides an estimate of the 
ratio of the global roughness exponent and the dynami- 
cal exponent : Q/ z « 0.26. The fit is computed for times 
between time t m i n and time t max . Before time t m i n the 
crack speed is not constant. After time t max , the rough- 
ness has saturated because of the reach of the system 
size. 

The lower line is a fit of the roughness measured for a 
small window size (I = 0.175 mm). It appears that w(l, t) 
increases like a power law as a function of the crack prop- 
agation time t even for small window sizes. The slope of 
the fit is 0.14 significantly larger than zero and gives an 
estimate of the /3* exponent. This unconventional depen- 
dence on time is an illustration of the anomalous scaling 
and differ from a Family- Vicsek scaling where the rough- 
ness is expected to be time independent for small window 
sizes. 

These two regimes are in good agrement with the 
anomalous scaling proposed in Eq.(||). A similar behav- 
ior has been observed for all specimens. 
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A. Local roughness exponent 

The local roughness exponent £/ oc is determined using 
four methods: the variable bandwidth methods : root 
mean square and max-min difference p|,p8|, the power 
spectrum and the averaged wavelet coefficient method 
p7| . Local roughness exponents Q oc have been deter- 
mined on profiles located far from the notch, i.e., at long 
times. Results on specimen s60-f are used as illustra- 
tions. Complete results for all specimens are provided in 
Table | 

In the root mean square method, the roughness w 
over a window I is expected from Eq.(|^) to scale at long 
enough time as 



(4) 



From Figure [| the local roughness exponent is : £/ oc = 
0.84 in the case of specimen s60-f . 

The max-min method consists of the computa- 
tion of h max (r), which is defined as the difference 
between the maximum and the minimum heights 
h within this window, averaged over all possi- 
ble origins x of the window pfj]: h max (r) — 

(Max{h(r')} x<rl<x+r - Min{h(r')} x<r i<x+r) x - For a 
self-affine profile, h max is expected to scale as: 



,(r) 



r-Cloc 



(5) 



where r is the width of the window along x-axis. For 
specimen s60-f , we measured a local roughness exponent: 
Ooc = 0.89. 

The third method is a calculation of the power spec- 
trum, i.e., the Fourier transform of the autocorrelation 
function (h(x + Ax)h(x)). The power spectrum scales for 
a self-affine profile as [Em : 



(6) 



where k is the wave factor. In Figure ^ we show a log-log 
plot of S(k) versus k for specimen s60-l. S(k) decays 
with a power law fc -2 78 which is consistent with (i oc — 
0.89. 

The last method used in this study is the averaged 
wavelet coefficient method |^7j]. This method consists of 
the average of the wavelet transform of the profile over 
the translation factor b. The averaged wavelet coefficient 
W[/i](a) scales as 



W[h](a) ~ a^+Cio 



(7) 



where a is the scale factor. 

The estimates of the local roughness exponents ob- 
tained with these four methods for all specimens are given 
in Table [| As shown in Table |, the values of the lo- 
cal roughness exponent Q oc calculated by the root mean 
square method, the max-min method and the power spec- 
trum decrease with the system size L. Only values ob- 
tained from the wavelet analysis seem independent of the 



system size. In the following, we show that this devia- 
tion is due to measurement and analysis biases and can 
be corrected. 

The reliability of the determinat ion of self-affine ex- 
has been 



ponents has already been studied 27 



It 



shown that several artifacts may introduce systematic er- 
rors for the estimation of the local roughness exponent. 
Two types of biases have to be distinghished: those which 
happen during the geometric measurement of the object 
and those which are relative to the method of signal anal- 
ysis. 

In our study, profiles are recorded with a needle moving 
along crack surfaces. For a similar type of measurement 
P,p8| , it has been shown that shape and volume of the 
needle can induce a geometric filter. When the tip of 
the needle is a half-sphere, the needle follows hills more 
correctly than sharp holes. Subsequently the exactness 
of the measured height is function of the surroundings. 
It has been found that an increase of the radius of the 
needle tip induces an increase of the measured roughness 
exponent (see [pll for more details). 

In the case of biases relative to the analysis methods, 
it has been found that the accuracy of the different meth- 
ods is sensitive to two parameters : the size of the system 
(number of recorded points) and the roughness exponent. 
In our study, the system size strongly evolves from small 
to big specimens : 360 to 2050 points. In Ref. p^j, tests 
on synthetic profiles generated with self-affine exponent 
between 0.8-0.9 show that the three methods underes- 
timate the self-affine exponent. The underestimation is 
larger when the system size decreases. The root mean 
square method is the most sensitive to this size effect. 

Likely both biases exist in this study. In order to eval- 
uate simultanously the influence of both flaws on local 
roughness exponents, synthetic profiles are simulated, fil- 
tered and analysed. Self-affine profiles are simulated nu- 
merically with a Voss construction f29[| for four values 
of self-affine exponent : 0.80, 0.85, 0.90 and 0.95. For 
each exponent, 100 independent profiles are generated. 
The horizontal step between two consecutive points is 
So = 6.25 /mi corresponding to the lower cutoff, i.e., the 
mean thickness of cell walls. Magnification of self-affine 
profiles corresponds to that measured on experimental 
profiles. The filter is an under sampling with a sphere 
of radius R — 25 /im (i.e. size of the experimental nee- 
dle) every four steps (S = ASo = 25/xm). The step S 
corresponds to the experimental step Air. Output expo- 
nents are obtained with the four methods (rms, max-min, 
power spectrum and wavelet analysis) and are given in 
Table |l| for different system sizes. From Table H, the cor- 
rected values of the experimental Q oc are estimated and 
given in brackets in Table |[ Average of the corrected 
values of Q oc obtained from the different methods gives 
the local roughness exponents : 0.87±0.07 for spruce and 
0.88±0.07 for pine. Results are consistent with those ob- 
tained for brittle materials where C/oc ~ 0.85 PflpTl and 
support the conjecture of a universal local roughness ex- 
ponent. 
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Our results are different from those obtained by Eng0y 
et al |3(| who studied the roughness of brittle fractures 
for different woods. The authors found a local rough- 
ness exponent Q oc — 0.68 which is characteristic of a 
two-dimensional fracture. Several reasons might explain 
this difference. First, direction of propagation crack was 
perpendicular to fibers while in our setup propagation 
is parallel to fibers. Second moisture content of tested 
specimens was around 4% which is significantly lower 
than that measured in our study (12%). Low moisture 
content induces micro-cracking in the radial-longitudinal 
and tangential-longitudinal planes of wood due to drying 
shrinkage. This mechanism of micro-cracking does not 
appear in mode I fracture. Micro-cracks induce preferen- 
tial paths for the macro-crack which modify the scaling 
properties of fracture surfaces. Third the experimental 
procedure was strongly different in the study of Eng0y 
et al since fracture propagation was unstable contrary to 
the stable propagation in the present work. 



B. Global roughness and dynamical exponents 

As discussed above, the existence of an exponent /3* ^ 
(see Fig. 2) indicates that an anomalous roughening 
is taking place. To obtain an accurate description of the 
anomalous scaling, we follow Jll],|24| and define the scal- 
ing function g(u) as g(l/ t x l z ) — w(l, t)/ From Eq.(||), 
g{u) is expected to scale like 



9{u) 



t -(C-Ci.„) 



if h«1 
if u » 1 



(8) 



The scaling function g is computed by data collapses from 
each profile of a complete crack map, [i.e. the set of pro- 
files that describe a single fracture). In Figures | and 
H we present the data collapses of g(u) for all the maps 
obtained for the three specimen sizes (L = 60, 30 and 
11.25 mm) of both spruce and pine. 

Fig.||(a) is considered as a good example of these data 
collapses. The quality of the collapse is used for the de- 
termination of the dynamical exponent z. The global ex- 
ponent £ is obtained from the fit of the scaling function. 
The time evolution of the height fluctuations at small 
scales is shown by the nonconstant behavior for m« 1. 
This regime is fitted by a power law g(u) oc u~ °' 76 . Us- 
ing our previous estimate of the local roughness exponent 
(loc = 0.84 we obtain the magnitude of the global rough- 
ness exponent £ = 1.60. In the particular case of sample 
s60-l shown in Fig. |](a) the best collapse is observed for 
the dynamical exponent z = 5.9. Note that estimates of 
the exponents Q oc = 0.84, £ = 1.60 and z = 5.90 are very 
consistent with fits obtained from Fig. |^: /?* = 0.13 and 
Ql z = 0.27 for this sample. 

Data collapses of all maps are presented in Fig.|^ and 
Fig.^J and are in good agreement with a scaling func- 
tion like (||). For both wood species, global rough- 
ness exponent and dynamical exponent are reported in 



Tab.|. As shown in Tab.|, global roughness exponents £ 
are independent of the system size. Average values are 
C = 1.60 ± 0.10 for spruce and C = 1-35 ± 0.10 for pine. 



IV. IMPLICATIONS OF ANOMALOUS SCALING 

According to Eq.(^), the roughness is expected to sat- 
urate only at times t 3> L z , i.e., when the correlation 
length ~ t x l z has reached the boundary length, 
imax oc L. In this regime the roughness magnitude scales 
with the system size for any window length even much 
smaller than the system size L: 



«;(/,*> L z ) ~ L^ 1 "" 



(9) 



We checked the linear relationship between t; max and the 
system size L by measuring £ max . From the evolution 
of the roughness w(l,t) with time (see Fig. ||), the sat- 
uration time t sa t is estimated. The correlation length 

x/z 

(,max is obtained using the relation £, max oc t sat . Values 
of £, m ax for the different maps are reported in Table Q. 
In figure fj], ^ ma x is plotted versus L for both woods. A 
linear relationship between ^ max and the system size L 
exists except in the case of spruce for the biggest sample 
size where the saturation regime is not clearly reached. 

In Figure [§ the ratio (w(l, t > ^ ax ) is plot- 

ted versus L for profiles at times t > (£, m ax) z for pine 
specimens. A power law L^~^ loc with exponents de- 
termined previously £ = 1.35 and Q oc = 0.80 is very 
consistent with data. It confirms the increase of the 
roughness magnitude with the system size L even for 
windows smaller than the system size. In Figure ^ 
(w(l <C £ max ? ))l versus L is also plot- 

ted (filled symbols) which is expected to be constant. 



V. CONCLUSIONS 

In this study we have shown that fracture surfaces of 
an anisotropic material like wood display an anomalous 
dynamic scaling of the crack roughness. From different 
specimen sizes, we have studied the size effects on rough- 
ness exponents. It appears that the global roughness ex- 
ponent is independent of the system size and different for 
both studied woods. We have obtained ( — 1.60±0.10 for 
spruce and £ = 1.35 ± 0.10 for pine. The local roughness 
exponent Q oc shows a deviation according to the system 
size. However, we argue that this deviation is due to a bi- 
ased estimate resulting from two independent effects: the 
number of sampled points and the local filtering result- 
ing from the needle shape during the roughness measure- 
ment. Errors due to these biases have to be considered 
and the corrected values of the local roughness exponents 
are 0.87±0.07 for spruce and 0.88±0.07 for pine. These 
results support the conjecture of a universal local rough- 
ness exponent for brittle materials. Moreover, we have 
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shown that it exists a linear relationship between the sys- 
tem size and the maximum correlation length ^ max . This 
relation induces a system size dependence in the rough- 
ness magnitude at saturation. 

Our results can be compared with a recent experiment 
in granite ]19| ] in which the exponents Q oc — 0.79 and 
£ = 1.2 were obtained. We suggest that the global rough- 
ness exponent, which seems to be dependent of material, 
may be a good candidate as an index for characterizing 
material properties. On the contrary, the local roughness 
exponent does not seem to change for different materi- 
als and might be universal. Up to our knowledge, the 
existing models of cracks are based on the assumption 
that cracks are truly self-affine,ie. C = Qoc- It is a ma- 
jor interest to find theoretical models of crack interfaces 
that could incorporate anomalous kinetic roughening in 
a simple way. 
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FIG. 1. Modified tapered double cantilever beam (TDCB) 
specimen subjected to mode I crack propagation. The crack 
plane is perpendicular to the tensile axis which corresponds to 
radial-longitudinal plane of wood (longitudinal direction be- 
ing the direction of crack propagation). Dimensions are given 
in mm. 

FIG. 2. Roughness (rms) w(l,t) vs. time for a spruce spec- 
imen 60 mm wide (s60-l) calculated over windows I of size 
ranging from I = 0.175 mm to I — 13.975 mm with size step 
Al = 0.100 mm. The continuous line (a) corresponds to the 
fit of the roughness between t m i n and t max , w(l,t) ~ t 13 * for a 
small window size / = 0.175 mm. /3„ = 0.14 is obtained. The 
continuous line (b) is the fit of data for large window size / = 
13.975 mm. Its slope 0.26 corresponds to f/a. 

FIG. 3. Roughness (rms) w(l) vs. I for the profile at the 
saturated time t 3> ^ m ax Z (specimen s60-l). The straight 
line corresponds to the power law w(l) ~ l^ 1 "" and gives a 
determination of the local roughness exponent £i oc = 0.84. 

FIG. 4. Power spectrum at time t 3at in the case of s60-l 
specimen. The straight line has slope -2.78 which is consistent 
with a power law fc - ( 2 ^°c +1 ) and a local roughness exponent 
0oc = 0.89. 
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FIG. 5. Data collapses for spruce specimens data in three 
different system sizes. Panel (a), (b), (c) display the data 
collapses of s60-l, s30-l and sll-1 specimens which have re- 
spectively a size L = 60, 30 and 11.25 mm. The non-constant 
behavior (i.e. non zero slope) at small values of l/t x ' z displays 
the dependence on time of the roughness magnitude, Scalings 
are in good agreement with the scaling function (Eq.(^)). 

FIG. 6. Data collapses for pine specimens of three different 
system sizes. Panel (a), (b), (c) display the data collapses of 
p60-l, p30-2 and pll-2 specimens having respectively a size 
L = 60, 30 and 11.25 mm. 



FIG. 7. Maximum self-affme correlation lengths ^max vs. 
the system size L for spruce (circle) and pine (square) spec- 
imens. Both spruce and pine show a linear relationship 
(dashed line) between £ m ax and L. The determination of 
£,max in the case of spruce and large system size is underesti- 
mated owing to the to little duration of the roughness map. 
The saturation transition was not clearly observable for this 
sample. 

FIG. 8. Size effect on the amplitude of the roughness over 
profiles at saturation, i.e., at times t > £, m ax z , for pine speci- 
mens. Upper symbols correspond to (w(l <C £max,t)/ l c,loa ) l 
versus L : Circle correspond to the specimen of L = 60 
mm, squares L = 30 mm, diamonds L — 22.5 mm and 
triangles L = 11.25 mm. Filled symbols are obtained for 

(W(l < £ m ax,t)/(l Cl °-= £max''~ Cl ° c )) l VerSUS L. 
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TABLE I. Description of the analysed specimens. Local roughness exponents are calculated using the root mean square 
method, the max-min difference method, the power spectrum, and the averaged wavelet coefficient analysis. Values in brackets 
are corrected from errors due to measurement and analysis biases. 





Specimen 


L 


Nb. of 


root mean 


Power 


max-min 


wavelet 


C 


z 


^max 


species 


label 


(mm) 


profiles 


square 


spectrum 




analysis 






(mm) 




s60-l 


60 


49 


0.84 (0.95) 


0.89 (0.86) 


0.89 (0.90) 


f.00 (0.96) 


1.60 


5.90 


3.90 




s60-2 


60 


49 


0.8f (0.88) 


0.85 (0.82) 


0.89 (0.90) 


0.9f (0.87) 


1.55 


2.40 


4.30 




s60-3 


60 


43 


0.84 (0.95) 


0.95 (0.93) 


0.87 (0.87) 


0.99 (0.95) 


1.60 


5.60 


not sat. 


spruce 


s30-f 


30 


48 


0.78 (0.84) 


0.83 (0.80) 


0.83 (0.79) 


0.92 (0.87) 


1.55 


3.50 


4.10 




s30-2 


30 


47 


0.79 (0.85) 


0.84 (0.81) 


0.82 (0.79) 


0.89 (0.84) 


1.60 


2.00 


4.20 




sff-f 


ff.25 


22 


0.73 (0.80) 


0.84 (0.85) 


0.83 (0.84) 


0.9f (0.84) 


1.55 


2.50 


1.35 




sff-2 


ff.25 


25 


0.77 (0.85) 


0.84 (0.85) 


0.83 (0.84) 


0.98 (0.91) 


1.60 


2.60 


1.45 




sff-3 


ff.25 


26 


0.79 (0.88) 


0.88 (0.93) 


0.84 (0.85) 


0.95 (0.88) 


1.55 


2.60 


1.40 


spruce 




0.88 ±0.05 


0.86 ±0.06 


0.85 ±0.07 


0.89 ±0.09 


1.60 ±0.10 








p60-f 


60 


45 


0.84 (0.95) 


0.91 (0.88) 


0.90 (0.93) 


0.97 (0.93) 


1.30 


1.90 


7.30 




p60-2 


60 


46 


0.8f (0.88) 


0.86 (0.83) 


0.88 (0.88) 


0.90 (0.86) 


1.35 


2.30 


not sat. 




p30-f 


30 


2f 


0.84 (0.95) 


0.87 (0.85) 


0.88 (0.88) 


0.99 (0.95) 


1.35 


2.20 


3.85 




p30-2 


30 


30 


0.80 (0.87) 


0.8f (0.79) 


0.87 (0.87) 


0.99 (0.95) 


1.30 


4.30 


not sat. 


pine 


p30-3 


30 


3f 


0.85 (0.95) 


U.83 (0.8UJ 


0.90 (0.95J 


0.9( (0.93J 


1.4U 


2. DO 


o on 
3.80 




p30-4 


30 


3f 


0.83 (0.92) 


0.88 (0.88) 


0.89 (0.90) 


1.01 (0.96) 


1.35 


3.90 


not sat. 




pff-f 


ff.25 


26 


0.75 (0.83) 


0.86 (0.88) 


0.83 (0.83) 


1.03 (0.96) 


1.35 


3.20 


1.40 




pff-2 


ff.25 


27 


0.75 (0.83) 


0.86 (0.89) 


0.83 (0.83) 


0.98 (0.91) 


1.40 


1.80 


1.80 




pll-3 


ff.25 


28 


0.75 (0.83) 


0.82 (0.79) 


0.84 (0.84) 


0.97 (0.90) 


1.30 


2.30 


not sat. 


pine 


p22-f 


22.50 




0.8f (0.89) 


0.84 (0.83) 


0.85 (0.83) 


0.99 (0.94) 






3.00 




p22-2 


22.50 




0.8f (0.89) 


0.86 (0.87) 


0.85 (0.83) 


0.94 (0.88) 






2.70 


pine 




0.89 ±0.05 


0.84±0.07 


0.87 ±0.05 


0.92 ±0.08 


1.35 ±0.10 







TABLE If. Tests of root mean square, power spectrum, max-min difference and wavelet analyses on undersampled and 
filtered synthetic self-affine profiles which model profiler recording (see text for details) . Four system sizes, in terms of number 
of points are considered. The accuracy of the exponents presented in this table is around 8%. 



System size 


256 pts 


512 pts 


1024 pts 


2048 pts 


self-afine exponent 


0.80 0.85 0.90 0.95 


0.80 0.85 0.90 0.95 


0.80 0.85 0.90 0.95 


0.80 0.85 0.90 0.95 


rms 


0.72 0.77 0.80 0.82 


0.74 0.78 0.81 0.83 


0.75 0.79 0.82 0.84 


0.75 0.79 0.82 0.84 


power spectrum 


0.82 0.84 0.87 0.89 


0.82 0.85 0.88 0.90 


0.83 0.87 0.89 0.91 


0.83 0.88 0.93 0.96 


max-min 


0.81 0.84 0.86 0.89 


0.83 0.85 0.88 0.90 


0.84 0.86 0.89 0.90 


0.84 0.86 0.89 0.91 


wavelet analysis 


0.87 0.92 0.97 1.02 


0.87 0.91 0.96 1.00 


0.85 0.90 0.94 0.99 


0.84 0.89 0.94 0.99 
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30 
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160 


320 
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560 





direction of crack propagation 



8.0 



6.0 I- 
4.0 - 

X 
CO 

E 

2.0 I- 



0.0 



5.0 



y 



y 



y 



y 



y 



y 



y 



y 



y 



y 



y 



y 



y 



y 



y 



y 



Us' 



y 



y 



y 



y 



/ □ 



y 



y 



y 



3* 



y 



y 



O 
O 



15.0 25.0 35.0 45.0 55.0 65.0 

L (mm) 



